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ABSTRACT 


This  report  contains  final  results  of  the  quantitative  characterizations  of  the  pore  and 
constituent  particle  morphologies  in  two  diffa^nt  variants  of  the  70S0-T74S1  six-inch- 
thick  plate  alloys.  The  purpose  of  those  characterizations  was  to  find  volume  firaction, 
size,  shape  and  spatial  distributions  of  both  phases  and  thdr  through-thickness  gradients. 
This  is  the  first  ever  characterization  of  this  scope  which  provides  unavailable  until  now 
data  necessary  for  the  predictions  of  the  size  distributions  of  the  &tigue  crack  initiating 
flaws,  and  for  the  development  and  test  of  all  microstructure  based  fatigue  life  prediction 
and  durability  models.  The  sizes,  sluq)es,  and  orientations  of  both  phases  were  analyzed 
using  the  linear  intooept  method.  Their  spatial  distributions  were  quantified  using  the 
nearest  neighbor  spacings  obtained  via  Voronoi  (Dirichlet)  tessellation  techniques.  Two 
types  of  intercept  lengths,  global  and  average  [5],  have  been  used  to  represent  sizes  of 
both  phases  on  differently  oriented  test  planes  and  at  different  locations  within  the  plates. 
The  distributions  of  those  size  measures  for  both  features  were  lognormal.  The  three 
dimensional  size  distributions  were  also  estimated.  They  were  obtained  from  the 
characteristics  of  the  computer  simulated  three  dimensional  model  microstructures  which 
had  the  same  two-dimensional  attributes  as  the  two  dimensional  sections  of  the  real  alloy. 
The  analysis  of  nearest  neighbor  spacing  indicated  that  both  the  pores  and  particles  have 
spatial  distributions  which  can  be  categorized  as  random  with  clusters.  In  addition,  an 
attempt  has  been  made  to  use  the  tessellation  cell  shapes  to  identify  types  of  spatial 
anisotropy  of  both  phases.  Results  are  currratly  used  in  the  development  of  the  extrone 
value  based  methodology  for  predicting  size  distributions  of  the  fatigue  crack  initiating 
flaws  and  in  the  deveiopmem  of  the  microstructure  based  stochastic  life  prediction  models. 
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INTRODUCTION 


The  most  important  microstructural  features  controlling  mechanical  behavior  of  the 
engineering  alloys,  and  aluminum  alloys  in  particular,  are  grains,  pores,  constituent 
particles  and  precipitates.  Since  the  sizes  and/or  spacings  of  those  features  are  usually 
non-unifcnm,  to  obtain  meaningful  corrdation  between  the  microstructure  and  mechanical 
properties  it  is  necessary  to  characterize  not  only  the  features'  average  sizes,  but  also  thdr 
size  and  spatial  distributions,  slu^)es,  orientations  and  volume  fractions.  The  purpose  of 
this  work  was  to  assemble  the  most  efficioit  quantitative  techniques  which  would  facilitate 
such  characterizations  and  to  test  them  by  quantifying  microstructural  gradients  in  the 
7050-T7451  A1  plate  alloy.  This  alloy  is  known  to  have  microstructural  gradients,  hence  it 
is  an  ideal  candidate  for  testing  gradient  characterization  methodologies.  The  grain 
structures  of  the  same  alloy  were  previously  studied  [1]  thus  the  focus  of  this  work  was  on 
the  characterization  of  the  large  second  phase  dispersions  such  as  pores  and  the 
constituent  particles. 

All  characterizations  of  the  second  phase  dispersions  require  two  types  of  measurements; 
(1)  the  measurements  needed  for  obtaining  average  individual  characteristics,  such  as  size, 
sh^  and  orientation,  and  (2)  the  measurements  ^lich  give  global  characteristics,  i.e. 
volume  fraction,  size  and  spatial  distributions  and  alignment.  The  sizes  of  the  second 
phase  dispersions  are  most  conveniently  measured  us'mg  the  linear  intercept  method  [1-4]. 
Obtained  size  distributions  can  be  then  analyzed  for  conformity  with  the  most  common 
standard  distributions,  such  as  normal  and  lognormal.  To  describe  shapes,  various 
arbitrarily  defined  shape  factors  can  be  used  [1-2,  5].  A  common  problem  associated  with 
these  factors  is  that  they  can  have  the  same  values  for  diametrically  different  shapes. 
However,  since  dispersions  in  engineering  materials  have  usually  similar  and  relatively 
simple  shapes,  their  shape  differences  can  be  adequatdy  described  using  the  shape  factors. 
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For  complex  shapes  several  more  precise  techniques,  such  as  Fourier  and  spectral  analysis, 
are  also  available  [6-8],  As  to  the  global  properties,  the  volume  fraction  can  be 
approximated  by  the  dispersions  area  fiaction  measured  on  plane  sections.  The  alignment 
can  be  characterized  using  one  of  many  measures  of  the  'degree  of  orientation'  exemplified 
by  the  Saltykov's  method  [9]  or  by  the  orientation  &ctor  proposed  earlier  by  the  authors 
[S],  The  measurements  of  the  spatial  distributions  of  the  second  phase  dispersions  will  be 
described  in  details  later  on. 

All  measurements  used  in  this  work  have  been  taken  firom  the  LS,  TL  and  ST  planar 
sections.  Fig.  1,  and  do  not  involve  conversion  fitctors.  These  particular  sections  were 
chosen  because  they  are  commonly  used  in  microstructural  characterization  of  the 
engineering  alloys.  To  accurately  and  completely  characterize  the  second  phase 
dispersions,  it  is  necessary  to  transform  their  two^imensional  (2-D)  measurements  to  the 
three-dimensional  (3-D)  size  and  spatial  distributions.  There  have  been  many  attempts  to 
develop  stereological  techniques  for  estimating  the  3-D  particle  size  distributions  from  the 
2-D  probes.  The  earliest  work  can  be  traced  back  to  1925  [10]  and  new  techniques  are 
proposed  each  year.  Underwood,  in  his  classic  monography  [2],  has  summarized  all  pre- 
1970s  methods.  Among  them  the  DeHofF's  approach  [1 1],  published  in  1962,  allowed 
for  obtaining  the  size  distribution  of  cylindrical  ellipsoids  from  the  measurements  made  on 
random  plane  sections.  Subsequently  DeHofif  and  Bousquet  [12]  developed  methodology 
for  obtaining  size  distribution  of  triaxial  ellipsoids  from  the  distribution  of  their  linear 
intercepts.  However,  the  application  of  both  methods  is  limited  because  they  assume  that 
all  spheroids  have  the  same  shape.  Other  stereolo^cal  models  published  before  1970  also 
assume  identical  particle  shapes  and  are  applicable  only  to  the  particles  with  the  most  basic 
shapes  [2],  Since  particles  in  the  engineering  materials  are  neither  simple-shaped  nor 
identical,  the  applications  of  the  above  models  to  that  class  of  materials  is  limited. 
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In  the  1980s,  a  number  of  stereological  methods  for  obtaining  3-D  distributions  for 
arbitrary  oriented  particles  with  different  shapes  have  been  proposed  [13-18].  The  point- 
sampled  intercept  method  (PSIM)  is  one  of  them  [19-21],  Its  advantage  is  that,  according 
to  Gundersen  et  a/.[13],  it  does  not  require  assumption  about  particle  shapes.  The  basic 
formula  used  in  PSIM  is; 

(1) 

which  relates  Vu,  the  volume-weighted  average  volume  of  the  particles  to  the  point- 
sampled  intercept  length.  Although  Eq.  (1)  ^ves  an  unbiased  estimate  of  the  average 
particle  volume,  it  is  important  to  emphasize  that  Ui.  is  an  unconventional  volume 
weighted  average  rather  than  the  number  average,  Un,  commonly  used.  Vv  is  calculated 
from  expression; 


n 


1=1 


(2) 


where  Vi's  are  volumes  of  individual  particles.  The  PSIM  method  requires  no  assumption 
about  the  particle  shapes,  but  it  does  require,  as  pointed  out  by  Gokhale  and  Vander 
Voort  [20],  that  the  structure  is  isotropic.  Moreover,  the  method  gives  only  the  average 
value  and  not  the  distribution. 


Sterio  [14]  proposed  another  technique  for  3-D  size  distribution  analysis.  In  his  disector 
method  a  counting  system  consisting  of  two  parallel  planes  separated  by  a  known  distance 
h  is  used.  An  unbiased  estimate  of  the  number  of  particles  in  the  specimen  can  then  be 
obtained  by  counting  the  number  of  particles  whose  transects  are  completely  or  partially 
inside  the  counting  system.  Combined  with  standard  stereological  relations  the  method 
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yields  the  mean  particle  sizes  on  the  disector  sections.  Three  dimensional  size  distribution 
can  be  estimated  using  the  similar  principle  from  the  sets  of  serial  sections.  This 
estimation  procedure  is  the  most  efficient  when  h  is  about  a  one-third  to  one-fourth  of  the 
mean  particle  height.  The  method  supposedly  ^ves  better  estimation  of  the  size 
distribution  than  PSIM.  However,  the  obvious  disadvantage  of  this  approach  is  that 
finding  the  particle  size  distribution  is  extremely  work  intensive.  For  small  particles  (in  the 
order  of  a  few  pm)  and/or  volume  firactions  less  than  1%,  the  measurements  are  also  very 
difficult.  For  example,  to  obtain  a  statistically  meaningful  value  of  the  average  pore  size  in 
the  70S0  alloy,  it  requires,  for  conventional  methods,  to  measure  pores  on  about  200 
fields.  In  the  case  of  the  disector  method  it  would  be  necessary  to  double  that  number.  In 
addition,  obtaining  required  two  parallel  plane  sections  few  pm  apart  is  extremely  difficult. 


A  very  general  method  for  estimating  3-D  distributions  of  sizes  and  shapes  of  creep 
cavities  fi’om  size  measurements  made  on  polished  plane  sections  has  been  recently 
proposed  by  T.  -S.  Liu  ef  a/.  [22].  The  problem  was  cast  in  the  form  of  a  double  integral 
equation  whose  solution  yields  the  size  distribution.  In  this  method  both  the  aspect  ratios 
and  the  orientations  of  spheroids  can  be  varied  and  the  method  reduces  to  the  DeHoff  s 
approach  [1 1]  for  particles  with  identical  shapes.  However,  the  shapes  of  the  particles  can 
be  either  oblate  or  prolate  but  not  both  at  the  same  time.  For  particles  lacking  rotational 
symmetry,  such  as  the  pores  and  constituent  particles  in  the  7050  alloy,  this  technique  is  of 
little  value.  Karlsson  and  Sandstrom  [23]  proposed  similar  method  for  evaluating  3-D  size 
distribution  of  inclusions  based  on  measurements  from  2-D  sections.  They  assumed  that 
inclusions  are  ellipsoidal  with  non-cylindiical  symmetry.  The  main  drawback  of  their 
approach  are  the  assumptions  that  all  inclusions  have  the  same  shape  factors  and  that  they 
are  all  aligned. 
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To  fully  characterize  second  phase  morphology,  the  size  measurements  have  to  be 
supplemented  with  the  characterizations  of  spatial  distributions.  The  spatial  distributions 
are  usually  classified  as  r^lar,  random,  clusto’ed  and  random  with  clusters.  The  regular 
distribution  of  points  on  a  plane  is  easy  to  recognize  as  the  points  are  arranged  in  a  distinct 
pattern.  To  characterize  such  distribution  it  is  necessary  to  specify  only  the  point  density 
and  to  identify  the  polygon  comprising  the  pattern.  The  random  distributions  of  points  can 
be  adequately  characterized  by  the  particle  area  density  and  the  parameters  of  the  Poisson 
distribution  describing  interparticle  spacings. 

Clusters  of  particles  are  more  difficult  to  describe.  The  advantages  and  disadvantages  of 
the  methods  developed  for  this  purpose  have  been  recently  reviewed  by  Vander  Voort 
[24].  Among  them  the  number  density  [25],  nearest  neighbor  spacing  [26-28],  Voronoi 
(Dirichlet)  tessellation  [25,  29-32]  and  image  amendment  methods  [33,24]  are  the  most 
often  used.  In  the  number  density  approach,  the  number  of  particles  per  unit  area,  is 
counted  in  successive  locations,  or  fields,  on  the  specimen  section.  For  each  field,  the 
mean  and  its  standard  deviation  are  found.  An  increase  in  the  standard  deviation  for  a 
given  field  indicates  a  greater  degree  of  clustering.  It  is  recommended  that  for  the  best 
result  the  size  of  the  field  should  be  equal  to  the  size  of  a  cluster.  In  the  image  amendment 
method,  all  particles  separated  by  a  distance  less  than  or  equal  to  a  chosen  critical  spacing 
are  fused  together.  The  number  of  fused  a^egates  represents  then  a  measure  of 
clustering.  The  area  fi’action  of  particles  in  each  cluster  and  the  size  of  the  cluster  can  also 
be  determined.  The  key  assumption  in  this  method  is  the  critical  interparticle  spacing  and 
the  method  accuracy  increases  with  the  increase  of  the  test  field  size.  In  the  nearest 
neighbor  method,  the  ratio  of  the  mean  nearest-neighbor  distances  for  the  real  and  the 
equivalent  random  microstructures  and  the  ratio  of  their  variances  are  calculated.  These 
ratios  are  then  used  to  determine  the  particle  spatial  distribution  type  [26].  In  the  Voronoi 
tessellation  method  the  distribution  of  the  tessellation  cell  areas,  and  the  cells'  area 
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variances  are  used  to  identify  the  particle  spatial  arrangement.  The  last  two  methods  are 
related  because  near-neighbors  can  be  uniquely  identify  after  constructing  Voronoi 
tessellation.  Both  methods  have  been  used  in  this  study  and  will  be  discussed  further  in  the 
sections  to  follow.  We  also  note  that  there  have  been  attempts  to  use  3-D  tessellation  to 
characterize  particle  spatial  distribution.  In  this  case  the  cells  form  space-filling  assemble 
of  irregular  polyhydra  and  the  algorithms  for  their  computation  are  readily  available  [34- 
36].  The  applications  of  those  3-D  tessellation  techniques  to  the  characterization  of  the  3- 
D  particle  spatial  distributions  have  been  studied  by  Saxl  [37],  Stoyan  [38],  Serra  [39]  and 
Santalo  [40]. 

The  70S0-T7451  plate  alloy  is  an  ideal  material  to  test  predictive  capability  of  all  described 
above  techniques.  This  alloy  has  been  recently  intensively  studied  with  the  focus  on  the 
improvement  of  its  fatigue  durability  [41-43].  Results  showed  that  although  most  of  the 
mechanical  properties  of  slightly  different  alloy  variants  are  practically  the  same,  there  are 
differences  in  the  fatigue  cracks  initiation  sites  in  different  variants  and  this  translates  into 
significant  differences  in  the  fatigue  lives  [44,45].  This  suggests  that  there  are 
microstructural  differences  between  variants  and  that  they  are  responsible  for  the 
differences  in  fatigue  durability.  In  addition  it  is  also  expected  that  each  alloy  variant  has 
different  through-thickness  microstructural  gradients.  In  this  work  the  differences  in  the 
7050-T74S1  plate  alloys  constituent  particle  and  pore  size  and  spatial  distributions,  which 
are  most  likely  responsible  for  fatigue  behavior,  have  been  quantify  using  described  above 
methodologies.  First,  the  linear  intercept  method  was  used  to  measure  the  average  pore 
and  particle  sizes  and  their  size  variations  at  different  locations  within  the  plates.  All 
resulting  2-D  size  distributions  were  tested  for  conformity  with  the  normal  and  lognormal 
distribution  types.  The  3-D  size  distributions  were  then  obtained  by  finding  a  3-D 
computer  simulated  model  microstructure  having  the  same  2-D  characteristics  as  the  true 
microstructures.  Secondly,  the  particle  and  pore  volume  fractions  were  estimated  by 
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measuring  their  area  fractions  on  plane  sections.  Thirdly,  the  shapes  and  oriratations  of 
the  particles  and  the  pores  were  characterized  using  the  methodology  proposed  eariier  by 
the  authors  [1,S].  Finally,  spatial  distributions  of  both  phases  were  identified  and 
characterized  using  a  new  method  based  on  the  linear  intercept  scan  of  the  tessellation 
cells. 

EXPERIMENTAL 

The  materials  investigated  in  the  study  were  two  variants  of  the  7050-T74S 1  plate  alloys 
designated  as  old  pedigree  ("high"  porosity)  and  new  pedigree  ("low"  porosity),  both 
obtained  from  Alcoa  Laboratories.  The  new  alloy  shows  improved  fatigue  life  as  a  result 
of  the  microstructural  refinements  obtained  through  proprietary  modifications  of  the 
composition  and  processing  conditions.  Samples  for  metallographic  examinations  of  both 
alloy  variants  were  taken  from  the  longitudinal  (TL),  transverse  (ST)  and  short  transverse 
(LS)  planes,  and  at  three  locations,  i.e.  at  the  surface,  quarter  depth  and  the  center  of  the 
as-received  plates.  The  locations  of  all  test  planes  are  shown  in  Fig.  1 . 

The  samples  used  in  metallographic  examinations  were  prepared  by  mechanical  grinding  to 
a  surface  finish  of 600  i^t  uring  wet  SiC  paper,  followed  by  polishing  on  a  felt  wheel  with 
1  pm  and  then  0.25  pm  diamond  pastes  nuxed  vdth  ethanol.  This  procedure  prodded  a 
clean  surface  with  little  particle  pull-outs.  Measurements  of  the  porosity  and  the 
constituent  particles  were  made  on  as-polished  specimens  to  avoid  alteration  of  particle 
shapes  and  number-densities  by  etching.  All  measurements  were  made  from  optical 
photographs  taken  at  magnification  900X.  A  few  SEM  micrographs  were  also  taken  and 
used  in  the  study  of  the  particle  shapes. 

The  magnification  900X  used  in  all  measurements  was  chosen  as  a  compromise  between 
efficiency  and  accuracy.  It  allowed  for  detection  of  a  total  of  3908  constituent  particles  in 
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I 

both  studied  alloys.  This  corresponds  to  a  total  scan  area  of 22.42  mm^  and  to  21 7 
I  particles  per  each  test  plane  shown  in  Fig.  1 .  The  number  of  detected  pores  was  3261, 

which  is  approximately  181  per  test  plane.  The  total  area  searched  for  the  pores  was 
I  44.4S  nun^  This  area  is  twice  the  area  scanned  for  the  constituent  particles  because  of  the 

pores'  lower  volume  fractions. 

All  micrographs  obtained  from  metallographic  examinations  were  digitized  and  saved  as 
computer  graphics  files  for  subsequently  analysis  using  our  custom-made  image  analysis 
program.  The  program  analyzes  the  images  using  the  linear  intercept  method  and 
automatically  calculates  the  average  intercept  lengths,  shape  and  orientation  indices  and 
aspect  ratios  of  all  second  phase  features.  It  also  finds  size  distributions  using  appropriate 
probability  plot  [46],  constructs  tessellation  ceil  networks  and  analyzes  particle  spatial 
distributions. 

To  find  second  phase  volume  fractions  it  was  assumed  that  they  were  equal  to  the  area 
fractions.  For  each  depth  location  the  volume  fraction  was  calculated  as  the  average  of 
the  area  fractions  on  TL,  ST  and  LS  planes.  For  each  plane  these  fractions  were 
calculated  as  the  ratio  of  the  number  of  pixels  inside  particles  to  the  total  number  of  pixels 
in  the  entire  field.  Alternatively,  the  volume  fractions  were  also  estimated  using  the  point- 
sampled  intercept  length  method  [19-21].  In  this  method,  the  volume  fraction  is  assumed 
equal  to  the  ratio  of  the  random  sampling  points  hitting  the  particles  to  the  total  number  of 
sampling  points.  However,  the  results  obtained  using  point-sampled  method  changed 
slightly  with  the  number  of  sampling  points,  therefore  this  method  was  abandoned. 

A.  Two-Dimensional  Size  Distributions 

The  sizes  of  the  pores  and  particles  on  the  plane  sections  were  measured  using  the  linear 
intercept  method.  In  this  method  the  average  intercept  length,  L^a),  of  a  two  dimensional 
object  at  a  scan  angle  a  can  be  calculated  as  [47] 
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///a) 


(3) 


where  A  is  the  object  area  and  ff/a)  is  the  tangent  height  defined  as  the  projection  of  the 
measured  object  on  the  direction  perpendicular  to  the  scan  lines.  The  average  intercept 
length  over  all  scan  angles,  L,p  is  defined  in  this  case  as; 

L,i  =  -fw2(a)L2(a)da  (4) 

^  0 

where  yt'/ocj  is  the  weight  factor  equal  to; 

W2(a)=  (5) 

ijH2(a)da 

’'o 

Jl 

For  convex  shapes,  j  H(a)da  is  equal  to  Lp,  the  perimeter  of  the  measured  object.  Thus 
0 

in  this  case 

(6) 

which  is  a  standard  equation  used  in  the  literature  [47].  To  study  the  directional  variation 
of  the  average  intercept  length,  we  also  define  a  normalized  intercept  length,  d(a),  as 

d(a)  =  d  =  |s^  (7) 

where  is  the  maximum  average  intercept  length.  The  d  vs.  a  plot  can  be 
subsequently  used  to  compare  object  shapes  [S]. 

To  compare  differences  between  the  pore  and  particle  sizes  of  the  new  and  the  old 
pedigrees  7050  alloys  the  standard  one  factor  ANOVA  (analysis  of  variance)  has  been 
performed.  In  this  analysis,  the  actual  F  value  is  compared  with  the  critical  F  value  for  5% 
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significance  level.  If  the  actual  F  value  is  smaller  than  the  critical  one,  then  the  null 
hypothesis  that  the  pore  or  particle  sizes  in  both  alloys  have  the  same  mean  could  not  be 
rejected.  This  means  that  the  average  pore  or  particle  sizes  of  the  two  pedigrees  are 
statistically  the  same  at  95%  probability. 


B.  Three-Dimensional  Size  Distributions 


The  transformation  of  2-D  measurements  to  the  3-D  size  distributions  has  been  carried  out 
assuming  that  all  second  phase  features  are  ellipsoids  with  different  sizes  but  identical 
aspect  ratios.  To  find  their  size  distributions  a  numerical  approacu  was  used.  A  computer 
program  developed  for  this  purpose  generated  a  model  microstructure  consisting  of  a  3-D 
spatial  distribution  of  ellipsoids  with  the  same  aspect  ratios  but  different  sizes.  Then  the 
program  sectioned  the  model  microstructure  along  TL,  LS  and  ST  planes,  estimated  the 
intercept  lengths  for  each  plane  and  compared  results  with  those  obtained  for  7050  alloy. 
When  the  differences  between  the  means  and  standard  deviations  of  the  2-D  size 
distributions  for  the  model  and  for  the  studied  alloy  were  minimum,  then  the  ellipsoids  in 
the  model  were  assumed  to  be  the  best  approximation  to  the  particles  in  the  real  alloy. 

The  measurement  of  the  goodness  of  fit,  R,,  used  by  the  program  had  the  following  form: 


l^i//.exp,i 


^  Lit  lexp.i 


?) 


(8) 


where  and  AOu,  are  the  differences  between  means  and  standard  deviations  of  L„  for 
the  experimental  and  simulated  microstructures  and  au,^^are  the  values  of 

mean  and  <;tandard  deviation  obtained  expenmentally.  The  data  for  ST,  LS  and  TL  planes 
correspond  to  i  =  1,  2  and  3  respectively.  The  definition  of  is  not  unique  and  other 
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forms  including  for  instance  weight  functions  added  to  both  terms  in  Eq.  (8)  could  also  be 
used. 

The  input  variables  to  the  3-D  model  were  (1)  dlipsoids  size  distribution  parameters,  (2) 
their  alignment,  and  (3)  their  spatial  distributions.  The  radii  of  ellipsoids  were  designated 
asa,  b  and  c  respectively  with  a  bdng  the  largest.  The  ratios  a/b  and  a/c  served  as  input 
parameters  and  they  could  vary  from  1  to  oo  to  accommodate  shapes  from  a  sphere  (a  =  6 
=  c)  to  an  arbitrary  ellipsoid.  The  size  of  each  ellipsoid  was  characterized  by  its  largest 
radius,  a,  and  it  was  assumed  that  those  sizes  obeyed  either  normal  or  lognormal 
distribution  The  maximum  radius  of  each  ellipsoid,  a  or  a[i],  was  obtained  from  the 
normal  standard  variable,  Z,  and  the  assumed  mean,  and  standard  deviation,  a, ,  using 
expression; 

aIil  =  p„+Z[i]a,  (9) 

In  case  of  the  lognormal  distribution,  all  size  dimensions  were  converted  into  their 
logarithmic  values.  Z[i]  was  calculated  using  the  method  by  D'Agostino  [48]  described  in 
the  Appendix. 

To  generate  a  random  dispersion  of  ellipsoids  in  a  given  volume,  the  coordinates  of  the 
center  of  each  ellipsoid,  (x,  y,  z),  were  independently  and  randomly  selected  within  the 
dimension  of  a  given  volume.  Each  location  was  then  checked  to  assure  no  overlapping 
between  ellipsoids.  If  overlapping  took  place  the  location  was  rejected  and  a  new  random 
location  generated  again.  As  a  result  the  spatial  distribution  of  the  ellipsoids  was  not 
exactly  random  but  pseudo-random  or  non-intersecting  random.  The  orientation  of  each 
ellipsoid  was  defined  by  three  Euler  angles  which  were  also  independently  and  randomly 
generated.  The  orientations  of  ellipsoids  in  the  model  could  be  varied  from  totally 
oriented  to  random.  The  orientation  distribution  was  assumed  to  be  known  a  priori,  and 
treated  as  an  input  variable  in  the  program.  The  spatial  randomness  of  the  simulated 
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mio’ostructures  was  tested  by  ecamining  a  collection  of  1,000  identical  spheres.  The 
theoretical  average  intercept  length  Lj  is  in  this  case  [2]; 


_  4V  _  4r 

"  S  “  3 


(10) 


vvdiere  V,  S  and  r  are  the  volume,  surfece  area  and  radius  of  the  sphere  respectively.  The 
value  of  Lj  obtained  for  the  simulated  microstructure  had  relative  error  of  only  2.9%. 


C.  Shanes  and  Alignment 


To  quantify  shapes  of  second  phase  features  the  shape  index,  SI,  equal  to  the  area  under 
the  d  vs.  a  curve  has  been  used  [1].  This  index  can  be  expressed  as 


SI 


±lh<S!l 

CL  •  L 

0  ^2.m«x 


(11) 


where  is  the  angular  period  of  the  d<t  curve.  As  demonstrated  elsewhere  [S],  SI  is 

independent  of  the  measured  object  size  and  orioitation.  It  has  values  between  zero  and 
one,  and  can  be  viewed  as  an  indication  of  how  drcular  is  the  object.  For  a  perfect  circle 
the  shape  index,  £/,  is  equal  to  one  and  its  values  for  typical  two-dimensional  shapes  can 
be  found  in  Ref  [5]. 


For  a  group  of  objects,  two  different  normalized  loigths,  the  average  d  and  the  globed  d 
[1,S]  can  be  used.  The  global  d-a  cui  obtained  uang  standard  intercept  length 
measuronents,  shows  the  change  of  the  object  normalized  intercept  length  with  directions. 
The  average  d<L  plot,  obtained  by  'aligning'  all  the  objects  before  scanning,  represents  the 
pure  object  shapes  as  the  orientation  effect  is  eliminated.  To  illustrate  the  difference  let's 
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assume  that  two  identical  elongated  particles  are  rotated  with  respect  to  each  other.  Their 
individual  d-a  curves  will  have  a  phase  difference  and  the  average,  represented  by  the 
global  d<L  curve,  will  reflect  both  the  shapes  and  the  orientations  of  both  particles. 
Consequently  the  global  shape  index,  SIg,  which  is  the  area  under  the  global  d-a  curve, 
represents  true  particle  shapes  only  if  the  particles  are  aligned.  Since  in  general  they  are 
not,  SIg  increases  with  increasing  particle  misalignment  and  approaches  unity  for  a 
statistically  random  structure.  The  "true"  average  shape  index  for  all  particles  can  be 
obtained  by  constructing  the  average  d-a  curve  for  particles  artificially  aligned  along  axes 
corresponding  to  We  shall  call  the  shape  index  resulting  from  such  curve  Slg^. 
Since  SI^  does  not  contain  the  information  about  particle  orientation,  comparing  it  with 
SIg  allows  for  the  deconvolution  of  the  shape  and  orientation  components.  This  is  done 
through  an  orientation  factor,  O,  defined  as; 


n= 


l-SI^ 


(12) 


For  randomly  oriented  features,  SIg  is  equal  one  and  will  be  zero.  For  a  perfectly 
oriented  structure  SIg  is  equal  to  5/^  and  Q  is  one.  The  orientation  factor,  Q,  allows 
thus  to  quantify,  on  a  scale  from  zero  to  one,  the  degree  of  particle  orientation.  This 
factor  is  used  to  describe  particle  orientation  distribution  in  our  3-D  simulation  program. 


Another  parameter  commonly  used  to  describe  the  degree  of  alignment  is  the  ratio  of  the 
total  length  of  oriented  lines  in  the  microstructure  to  the  total  length  of  all  lines,  Q;  ^ 

This  ratio  was  introduced  by  Saltykov(6]  and  in  our  case  it  is  equal  to  [1]; 


O  = _ UcL=o)-Ly(a^9o) 

L/a=o)  +  0.57 1  L/a=9o) 


(13) 
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Clf  2  is  zero  for  a  random  structure  and  approaches  one  for  structure  with  highly  oriented 
and  elongated  particles.  However,  it  is  not  guaranteed  that  the  structure  is  random 
Qj  2  is  zero  because  eq.  (13)  is  strictly  valid  only  for  partially  oriented  lines.  Therefore,  we 
prefer  to  use  Q  in  Eq.  (12),  which  uses  the  average  values  over  all  the  scan  angles  as  a 
more  precise  measure  of  particle  alignment. 

D.  Spatial  Distributions 

The  spatial  distributions  of  the  second  phase  particles  have  been  studied  with  the  hdp  of 
the  Voronoi  tessellation  constructions[2S-32].  Two  dimensional  Voronoi  cells  are  plane¬ 
filling  polygons  vnth  the  sides  of  each  polygon  constructed  as  the  bisectors  of  the  lines 
connecting  centers  of  the  neighboring  particles.  This  way  each  polygon  contains  one 
second  phase  particle  and  the  ratio  of  the  polygon  area  to  its  particle  area  is  defined  as  the 
"local  area  fraction".  With  the  help  of  tessellations  the  spatial  distribution  of  particles  can 
be  analyzed  by  studying  the  particle  local  area  fractions,  average  cell  areas  and  their 
standard  deviation  and/or  average  numbers  of  cell  sides  and  their  standard  deviation,  etc. 
[25-32].  However,  the  obtained  results  are  approximations  as  ideally  three  dimensional 
Voronoi  tessellations  should  be  used  for  spatial  charactoizations.  Another  method  for 
classifying  the  spatial  distribution  of  the  second  phase  particles  is  the  nearest  neighbor 
spacing  method  [26]  which  utilizes  two  indexes,  Q  and  R,  defined  as  follows; 

Q  _  Observed  Nearest  Neighbor  Distance 
Expected  Nearest  Neighbor  Distance 


I 

I 


^  _  Observed  Nearest  Neighbor  Distance  Standard  Deviation 
Expected  Nearest  Neighbor  Distance  Standard  Deviation 


(15) 
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where  the  expected  values  are  calculated  by  assuming  a  random  structure  with  the  same 
point  density  as  the  experimental  one.  Note  the  above  definition  of  R  is  based  on  the 
standard  deviations  instead  of  the  variances  as  originally  proposed  by  Schwartz  and  Exner 
[26].  However,  the  two  definitions  converge  to  the  same  exact  critical  values,  R  ==  1, 
under  the  same  conditions.  To  identify  different  types  of  spatial  distribution  the  following 
guidelines  can  be  used:  (1)  for  the  clustered  distribution,  Q  <  1  and  R  <  1;  (2)  for  r^lar 
distribution  Q  >  1  and  R  >  1;  (3)  for  random  distribution  with  clusters  Q  <  1  and  R  >  1, 
and  (4)  for  totally  random  distribution  Q  =1  and  R  =  1 . 

In  addition  to  using  Q  and  R  to  identify  spatial  distribution  patterns,  a  tessellation  cell 
shape  method  has  been  tried.  In  this  method,  the  sizes  and  shapes  of  the  tessellation  cells 
are  analyzed  using  the  same  methods  as  those  for  the  second  phase  dispersions.  The 
advantage  of  this  approach,  originally  suggested  by  Wary  et  al.  [27],  is  that  it  allows  for 
the  detection  of  the  fluctuations  in  the  cell  orientations.  By  using  this  technique,  all 
methodologies  developed  for  characterizirg  2-D  objects  can  be  applied  to  characterize  the 
polygonal  tessellation  cells.  Consequently  for  a  totally  random  set  of  points  representing 
second  phase  centers,  the  resulted  Voronoi  tessellation  cells  should  also  be  randomly 
oriented  and  their  shape  index,  SI,  should  be  equal  to  one.  If  the  point  distribution  is 
regular,  for  example  rectangular  or  hexagonal,  the  shape  indexes  are  smaller  than  one  and 
their  values  are  also  known  [5].  The  d-a  curves  for  the  tessellation  cells  constructed  on 
particles  or  pores  can  be  thus  compared  with  those  for  the  random  and/or  the  regular 
patterns  to  identify  their  spatial  distribution  type.  Since  most  results  published  so  fiu* 
correlate  the  cell  area  distributions  to  the  point  distributions  [29,  32],  one  can  use 
tessellation  cell  shape  method  to  gain  additional  information  on  the  particle  spatial 
distribution.  The  tessellation  cell  shape  method  is  particularly  useful  for  analyzing  the 
orientations  of  the  second  phase  clusters  which  cannot  be  detected  by  the  other  two  spatial 
distribution  methods.  This  is  because  in  an  ideal  clustered  structure  any  location  and  /or 
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orientation  changes  of  a  whole  cluster  will  not  affect  the  values  of  Q  and  R  because  the 
nearest  neighbor  spacii^  remains  the  same.  The  tessellation  cell  shapes,  on  the  other 
hand,  are  sensitive  to  any  anisotropy  in  the  cluster  orientations  spacings  and  cells  d-a 
curves  and  SI^s  indexes  can  be  used  to  analyze  those  inhomegeneities.  As  a  result  the 
information  on  the  spatial  distributions  provided  by  the  cell  shape  method  and  the  Q-R 
index  method  are  complementary;  the  former  ^ves  information  about  orientation  isotropy 
and  the  latter  characterizes  the  spatial  distribution  type.  All  above  methods  have  been 
used  to  analyze  pore  and  the  constituent  particle  patterns  in  this  work. 

RESULTS  AND  DISCUSSIONS 

A.  Pores 

1.  Volume  fraction 

The  average  pore  volume  fractions,  the  number  of  detected  pores  and  the  number  of  fields 
used  in  the  measurements  are  listed  for  each  test  plane  in  Table  1 .  As  already  explained, 
the  volume  fractions  were  assumed  equal  to  area  fi^ions  and  their  average  values  were 
0. 135%  for  the  old  and  0. 103%  for  the  new  pedigree  alloys.  The  porosity  level  in  the  new 
7050  alloy  has  been  thus  reduced  by  24%.  The  through-thickness  porosity  volume 
flection  gradients  are  shown  in  Figure  2,  with  each  data  point  representing  the  average  for 
TL,  LS  and  ST  planes  at  a  given  plate  location.  As  expected,  the  average  pore  volume 
finctions  are  the  highest  at  the  center  and  the  lowest  near  the  surface  for  both  studied 
plates.  However,  the  porosity  levels  are  practically  the  same  for  both  pedigrees  at  the 
surface  and  quarter-depth  locations  but  are  higher  at  the  center  regions  for  the  old  alloy. 
The  reduction  of  the  average  porosity  level  in  the  new  alloy  comes  then  mainly  from  the 
center  region. 
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2.  Size,  shape  and  alignment 


A  typical  pore  and  shapes  of  201  pores  found  on  the  ST  plane  at  the  coiter  of  the  old 
plate  are  shown  in  Figs.  3  and  4  respectively.  The  results  of  quantitative  charact«rizations 
of  these  and  all  other  pore  sizes,  shapes  and  alignments  are  summarized  in  TaUe  1.  The 
avenge  pore  dzes  are  represented  in  Table  1  by  the  average  intercept  length, 
measured  at  scan  angles  from  0  to  180**  in  incranents  of  S**.  To  show  the  change  of  pore 
intercept  lengths  with  directions  the  curves  have  been  used.  Figs.  S  and  6.  These 
curves  are  ellipsoidal  and  they  represent  particle  profiles  on  different  plane  sections  in  3-D 
space.  The  information  provided  by  Ljj  and  Lj-a  plots  are  complimentary.  Ljj  represents 
the  overall  average  pore  size  dianges  while  Lf<x  plots  show  the  variation  of  the  pore 
sizes  with  orientation  and  location  within  the  plate.  The  average  values  of  Ljj  were  3.39 
pm  for  the  old  and  3.47  pm  for  the  new  alloy.  The  through-thickness  gradients  of  Ljj  are 
shown  in  Figure  7.  For  the  old  70S0  variant,  the  average  pore  size  is  the  smallest  near  the 
plate  surface  and  increases  toward  the  plate  center.  For  the  new  pedigree,  the  anallest 
pore  size  is  also  near  the  plate  suifoce  but  there  is  almost  no  size  change  between  quarter- 
depth  and  the  cento*.  Considering  the  variation  of  pore  volume  fhu^tions.  Fig.  1,  it 
appears  that  both  studied  alloys  have  strong  microstiuctural  gradients  and  the  plate 
centers  has  the  largest  pores. 

Cumulative  2-D  pore  size  distributions  for  both  investigated  alloys  have  been  compiled  in 
Figure  8.  The  pore  size  measure  used  in  the  figure  is  as  before  L/,,  the  average  intercept 
length.  The  plots  in  Fig.  8  show  both  the  size  ranges  and  the  differences  between  the  pore 
sizes  at  different  locations  for  each  pedigree.  Comparison  of  the  distributions  suggests 
that  the  pores  in  the  new  alloy  are  not  only  smaller  but  also  have  narrower  distribution 
than  those  in  the  old  one.  To  further  quantify  the  size  distributions  they  were  fitted  with 
both  the  normal  and  the  lognormal  distributions,  using  appropriate  probability  plots[46] 
and  ^ploying  the  corrdation  coefficient,  to  measure  the  goodness  of  fit.  In  ail  cases 
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the  pore  sizes  followed  the  lognormal  distribution  with  dose  to  one.  The  obtained 
distribution  parameto^  are  listed  in  Table  2  and  an  example  of  a  t^ical  size  frequency 
distribution  and  the  lognoimal  fit  are  shown  in  Figs.  9  and  10. 

The  3-D  pore  size  distributions  were  estimated  using  described  earlier  numerical  method 
and  assuming  that  pores  could  be  represented  by  triaxial  dlipsoids.  The  results  are 
presented  in  Figure  1 1  together  with  the  original  experimental  data.  The  3-D  size 
measures  used  in  Fig.  1 1  are  the  maximum  int^cept  lengths,  L/jj^ ,  in  the  plane 
containing  dlipsoid  center  and  the  nuyor  axes.  As  in  the  case  of  the  2-D  size  distribution, 
the  lognormal  distribution  gives  a  better  match  b^ween  the  simulated  and  experimental 
results  than  the  normal  one.  The  means,  standard  deviations  and  aspect  ratios  of  the  3-D 
lognormal  size  distributions  for  both  alloys  are  listed  in  Table  3.  As  expected.  Fig.  1 1,  the 
maximum  average  intercept  loigths  of  the  3-D  dlipsoids  are  larger  than  the  2-D 
experimental  results.  This  is  because  for  a  ^vm  ellipsoid  any  2-D  intercept  length  is 
always  snudler  than  or  equal  to  its  3-D  Lj,  ^.  Consequently  the  experimental  2-D  data 
should  iq)proach  the  3-D  values  for  larger  pore  sizes  as  it  is  iiKleed  the  case  in  Fig.  1 1 . 

The  goodness  of  fit  of  the  estimated  3-D  size  distribution  was  evaluated  using  the  value  of 
the  normalized  errors,  >^ch  are  also  listed  m  Table  3.  Altemativdy,  the  error  can  be 
evaluated  by  conqraring  the  2-D  experimental  size  distributions  with  their  equivalents  for 
the  simulated  microstiucture.  An  example  of  such  comparison  is  shown  in  Figure  12 
where  the  two  resulting  curves  are  very  close  over  the  entire  size  range.  The  fits  could  be 
further  improved  if  otho*  distributions,  for  example  Gumbel  or  three  parameter  lognormal 
and/or  dififo'ent  aspect  ratios  for  dlipsoids  were  used.  Woric  on  the  implementation  of 
these  improvements  is  in  progress.  We  also  realize  that  the  format  of  the  normalized  error 
R,  is  not  unique  and  that  the  results  will  be  diffo’ent  for  differently  defined  R^'s.  Despite 
these  limitations,  the  obtained  approximations  of  the  3-D  pore  sizes  have  sufficioit 
accuracy  for  use  in  the  mechanical  behavior  models. 
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To  compare  the  pore  sizes,  Ljj,  in  the  new  and  the  old  7050  alloys  the  one  &ctor  ANOVA 
has  been  performed.  The  experimental  F-value  is  in  this  case  equal  to  3. 177  and  it  is 
smaller  than  the  critical  F-value  of  3.844  for  tlw  S%  confidence  level.  Consequently,  the 
null  hypothesis  that  the  pore  sizes  for  both  pedigrees  have  the  same  mean  could  not  be 
rejected,  and  the  average  pore  suzes  of  the  two  pedigrees  are  statistically  the  same  at  95% 
probability.  However,  the  tails  of  the  size  distributions  are  very  different,  especially  at  the 
larger  size  regime.  For  example,  there  are  23  pores  in  the  old  pedigree  with  diameters 
larger  than  10  pm  but  only  9  for  the  new  pedigree.  Since  only  the  largest  pores  act  as  the 
fiitigue  crack  initiation  sites,  it  is  the  upper  tail  of  the  size  distribution,  not  the  average 
value,  that  controls  the  fiitigue  behavior.  This  illustrates  the  importance  of  the  complete 
quantitative  microstructural  analysis. 

The  pore  shapes  have  been  characterized  using  the  d-a  plots  [1,5],  Figs.  13  and  14.  Fig. 
13  shows  the  global  d  and  the  average  d  curves  in  rectangular  coordinates.  Fig.  14  shows 
the  change  of  average  d  on  polar  plot.  Note  that  Fig.  14  represents  the  average  outline 
profile  of  the  pores.  The  shape  index,  57,  calculated  from  d-a  plots  in  Fig.  13  for  the 
pores  in  the  center  region,  ST  plane  of  the  old  alloy,  was  0.82.  The  average  shape  indexes 
for  pores  on  all  test  planes  were  0.816  and  0.831  for  the  new  and  old  alloys  respectively. 
This  means  that  the  average  shape  of  the  pores  on  2-D  sections  is  ellipsoidal  with  the 
average  ratios  of  the  shortest  to  the  longest  axes  of  0.816  and  0.83 1  for  new  and  old 
alloys  respectively.  The  corresponding  average  pore  aspect  ratios  defined  as  the  ratios  of 
the  average  at  0  to  90°,  were  1.45  and  1.51.  Since  the  pores  were  irregular,  the  shape 
index  57,  which  contains  information  from  all  the  scan  angles,  was  a  better  descriptor  of 
the  shapes  than  the  aspect  ratio.  However,  the  results  fi-om  both  methods  suggest  that  the 
pore  morphology  in  both  studied  alloys  were  sinular  as  the  relative  difference  betweoi  57's 
was  less  than  2%  and  for  the  aspect  ratios  it  was  only  4%.  The  change  of  pore  shapes 
with  location  within  the  plate  is  shown  in  Figure  1 5.  It  is  evident  that  the  changes  are 


20 


small  and  that  the  pores  in  the  old  alloy  are  slightly  more  elon^ed  than  their  counterparts 
in  the  new  alloy.  We  attribute  these  differences  to  different  thermo-mechanical 
treatments. 

The  orientation  fectors,  Q,  for  the  pores  on  different  sections  are  listed  in  Table  1.  They 
vary  from  0. 1 1,  which  corresponds  to  almost  random  structure,  to  0.6S  which  indicates 
partial  alignment.  The  average  values  are  0.38  for  the  new  and  0.32  for  the  old  pedigrees. 
Since  Q  in  Eq.  (12)  increases  as  the  pores  become  more  and  more  circular  and/or  more 
randomly  oriented,  the  alignment  of  pores  in  the  studied  alloys  can  be  considered  as 
insignificant.  The  change  of  the  average  O  values,  for  ST,  TL  and  LS  planes  with 
location  are  shown  in  Figure  16.  Results  suggest  slightly  larger  degree  of  alignment, 
particularly  at  the  surface  region,  in  the  new  alloy.  These  average  Q  values  were  used  as 
input  data  in  the  numerical  estimation  of  the  true  ellipsoid  sizes. 

The  results  of  the  characterizations  of  the  pore  shapes  and  orientations  provide  important 
information  to  the  fatigue  life  prediction  models.  Since  the  largest  pores  act  as  the  fatigue 
crack  initiation  sites,  their  shapes  and  sizes  determine  the  maximum  dimension  of  the 
'initial  flaw*.  The  alignment  factor,  Q,  offers  additional  information  on  the  relative 
orientation  of  the  pore  axis  with  respect  to  the  axis  of  the  applied  stress.  Both  of  them  are 
important  factors  controlling  frtigue  lives. 

B.  Constituent  Particles 

1.  Volume  fraction 

The  7050-7451  plate  alloy  contains  two  different  types  of  particles:  Al7Cu2Fe  and  MgjSi. 
Since  both  types  were  found  at  the  fatigue  crack  initiation  sites  they  were  characterized 
without  distinction  to  their  chemical  compositions  using  the  same  methods  as  for  pores. 
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The  average  particle  volume  fractions,  the  number  of  measured  particles  and  the  numbo' 
of  viewing  fields  used  in  the  measurnnmts  are  listed  in  Table  4.  The  obtained  constituent 
particle  area  fi'actions  were  0.61%  for  the  new  and  0.70%  for  the  old  alloy.  This 
difference  is  the  most  significant  improvemoit  in  the  new  70S0-T4S 1  plate  alloys  which 
came  out  from  our  microstructural  analysis.  The  through-thickness  particle  volume 
fraction  gradients  have  been  plotted  in  Fig.  2.  They  indicate  that  new  pedigree  contains 
fewer  particles  at  both  center  and  quarter  depths  than  the  old  one  but  the  opposite  is  true 
at  the  surface.  Despite  that  the  new  alloy  has,  on  average,  fewer  particles  than  the  old 
one. 

2.  Size,  shape  and  alignment 

A  typical  collection  of  constituent  particles  found  on  the  ST  plane  at  the  center  of  the 
plate  in  the  old  alloy  is  shown  in  Fig.  17.  The  average  2-D  sizes  for  constituents, 
measured  on  ail  nine  sections,  were  4.02  pm  for  the  new  pedigree  and  4.77  pm  for  the  old 
alloy.  Table  4.  The  average  "size  profiles"  of  the  constituent  particles  are  shown  in  Figs. 

18  and  19.  The  apparent  particle  sizes  are  always  the  smallest  on  the  ST  planes  due  to  the 
elongated  particle  shapes  resulting  fi-om  rolling.  The  null  hypothesis  that  the  particle  sizes 
for  the  new  and  old  alloy  were  the  same  had  to  be  rejected  since  the  resulting  experimental 
F-value,  34.272,  is  much  larger  than  the  critical  F-value  of  3.844  for  5%  significance  level. 
Therefore  the  average  particle  size  in  the  old  alloy  is,  at  95%  probability  level,  larger  than 
that  in  the  new  one.  The  change  of  partide  sizes  with  depth  for  both  alloys  is  shown  in 
Fig.  7.  The  size  increases  toward  the  plate  center  and,  for  example,  the  average  sizes  for 
the  new  alloy  are  3.25,  3.87  and  4.93  pm  at  the  surface,  quarter  depth  and  the  center  of 
the  plate  respectively.  The  differences  are  most  likely  due  to  the  temperature  variation 
during  the  thermal-mechanical  processing.  The  cumulative  particle  2-D  size  distributions, 
typical  size  frequency  distribution  and  the  lognormal  fit  are  shown  in  Figs.  20,  21  and  22 
respectively.  As  for  the  pores,  the  particle  sizes  also  have  lognormal  distribution  with 
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close  to  one.  The  parameters  of  the  obtained  lognormal  fit  are  listed  in  Table  S  and 
characteristics  of  the  calculated  true  3-D  size  distribution  are  summarized  in  Table  3. 


The  d-a  plots  for  constituent  particles  are  shown  in  Fig.  23  and  24.  Their  analysis 
revealed  that  particles  in  the  old  alloy  had  average  shi^  index  SI  equal  0.74  and  aspect 
ratio  of  1.86.  For  the  new  alloy  these  values  were  0.68  and  2.28  respectively.  Table  4. 
This  means  that  the  constituents  are  elongated,  ellipsoidal  in  sh^  and  aligned  with  the 
rolling  direction[l].  However,  we  have  also  observed  instances  of  isolated,  round  and 
slightly  oval  constituent  particles  as  well.  The  through-thickness  shape  index  variations 
are  shown  in  Fig.  IS.  They  suggest  that  particles  in  the  old  alloy  are  more  circular 
particularly  at  quarter  depth  and  surface  regions.  The  average  Q  values  for  the 
constituent  particles.  Table  4,  are  0.S2  and  0.43  for  the  new  and  old  alloys  respectively 
indicating  a  moderate  alignment.  The  change  of  O  from  the  surface  to  the  center  of  the 
plate  is  shown  in  Fig.  16.  The  results  show  that  the  alignment  exists  only  at  the  sur&ce 
and  center  regions  of  the  plates.  Comparing  characteristics  of  the  pores  and  the  particles 
in  both  alloys  it  appears  that  the  particles  have  higher  volume  fraction  and  larger  sizes  and 
that  they  are  more  elongated  and  aligned  than  pores. 

C.  Spacing  and  Clustering  of  Pores  and  Particles 

The  particle  and  pore  spatial  distributions  have  been  quantified  using  the  tessellation 
method.  Figure  25  shows  an  example  of  the  Voronoi  tessellation  cell  construction  for  the 
pores  on  the  ST  plane  at  the  center  location  of  the  old  alloy.  Table  6  lists  some  of  the 
characteristics  obtained  from  the  analysis.  The  most  important  among  them  are  local 
volume  fi^ions  and  near  and  nearest  neighbor  distances  —  these  quantities  can  be 
obtained  unambiguously  only  from  the  tessellations.  The  observed  nearest  neighbor 
distance  in  this  case  was  48.2  pm  with  standard  deviation  of  33.2  pm.  The  expected 
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counteqiarts  for  the  random  structure  with  the  same  point  density  are  5 1 .86  and  27. 1 8  (tm 
respectively.  This  gives  Q  =  0.929  and  R  =  1.224  which  suggests  that  the  distribution  of 
the  pores  can  be  classified  as  random  with  clusters.  For  the  constituent  particles  at  the 
same  location  the  observed  nearest  neighbor  distance  was  2S.2  pm  with  standard  deviation 
19.2  pm  while  expected  values  are  27.64  and  14.45  pm  respectively.  This  gives  Q  = 

0.879  and  R  =  1 .281  which  again  suggests  that  the  distribution  is  random  with  clusters. 
Above  information  about  the  second  phase  spatial  distributions  on  all  sections  at  the 
centers  of  the  plates  have  been  compiled  in  Table  7.  In  all  cases  the  spatial  distributions 
are  random  with  clusters.  It  is  also  possible  to  carry  the  spatial  characterizations  further 
by  considering  the  second  phase  cluster  density,  the  number  of  the  features  per  cluster,  the 
cluster  size,  etc.  and  we  are  currently  working  on  the  implementation  of  methodologies 
allowing  for  such  characteristics.  We  again  emphasize  that  ideally  one  should  use  3-D 
tessellations  to  study  the  spatial  distribution.  As  pointed  out  earlier,  the  main  difficulty  in 
this  case  is  acquisition  of  the  data  on  the  3-D  particle  locations,  which  makes  the  method 
impractical. 

An  example  of  the  results  obtained  fi'om  the  shape  analysis  of  the  tessellation  cells  is 
shown  in  Fig.  26.  The  ceils  analyzed  are  those  from  Fig.  25  and  the  method  used  is  the 
same  as  the  one  utilized  for  characterizing  shapes  of  pores  and  particles.  Fig.  26  shows 
the  d-a  curve  obtained  for  the  cells  and,  for  comparison,  the  curves  for  the  random  and 
rectangular  point  patterns.  The  SIg  for  the  tessellation  is  0.963  while  it  is  1  for  the  random 
and  0.760  for  the  rectangle  patterns.  The  shape  index,  Sl^^,  is  in  this  case  equal  to  0.736 
which  gives  orientation  factor  Q  of  0. 138.  The  average  cell  aspect  ratio  is  1 .872.  As 
discussed  earlier,  a  high  value  of  5/^  combined  with  a  low  value  of  Q  indicates  near 
isotropic  cell  structure.  The  values  of  aspect  ratio  and  SI^  suggest  that  individual  cells 
are  in  this  case  elongated  and  ellipsoidal.  Work  on  the  refined  version  on  the  cluster 
characterization  methodology  is  in  progress. 
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CONCLUSIONS 


•  The  volume  fractions,  sizes,  shapes  and  orientations  of  pores  and  constituent  particles 
in  two  different  variants  of  70S0-T74S1  A1  plate  alloy  have  been  quantified.  The 
characterizations  included  measurements  of  not  only  the  average  values  of  all 
microstructural  descriptors,  but  also  their  distributions  for  TL,  ST  and  LS  planes  at  the 
surface,  quarter  depth  and  at  the  center  of  each  plate. 

•  The  lognormal  distribution  has  been  found  to  give  the  best  description  of  the  size 
distribution  for  both  pores  and  particles.  The  average  pore  sizes  were  statistically  the 
same  in  both  studied  alloys  but  the  average  particle  size  in  the  old  alloy  is  larger  than 
that  in  the  new  one.  The  new  alloy  has  improved  microstructures  which  contains  24% 
less  pores  and  has  smaller  particles  than  the  old  one. 

•  In  both  alloys  the  pores  are  elongated  with  the  average  aspect  ratio  of  1 .50.  For 
constituent  particles  the  average  aspect  ratio  is  2.0.  The  spatial  distributions  of  both 
pores  and  particles  were  characterized  using  the  nearest  neighbor  spacing  and  the 
Voronoi  tessellation  methods.  In  all  cases  they  can  be  classified  as  random  with 
clusters. 

•  A  numerical  method  for  estimating  the  3-D  size  distributions  from  the  experimental  2- 
D  measurements  made  on  plane  sections  has  been  developed.  The  method  assumes 
that  particle  sizes  and  shapes  can  be  modeled  as  collection  of  triaxial  ellipsoids  with 
random,  partially  aligned  or  totally  aligned  orientations.  The  method  has  been  used  to 
obtain  3-D  size  distributions  of  pores  and  constituent  particles  in  the  7050-T7451  A1 
plate  alloy. 
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•  The  shapes  of  tessellation  cells,  built  on  particle  or  pore  centers,  have  be^  analyzed  to 
quantify  the  second  phase  spatial  distributions.  The  results  have  been  compared  with 
those  for  random  and  various  regular  point  sets  to  provide  additional  information 
about  the  pore  and  particle  spatial  distribution  patterns. 
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APPENDIX 


The  standard  random  variable,  Z,  can  be  i^)proxiniated  as[46] 

Z  =  sigriiF^  (jc)  -  0.511.238/(1 + 0.0262/)] 

where 


(Al) 

(A2) 

(A3) 

(A4) 


where  N  is  the  sample  size  and  i  is  the  observation  number.  To  obtain  probability  plots 
one  needs  to  plot  the  standard  variable  Z  has  been  plotted  as  a  function  of  a  random 
variable  X;  for  the  normal  or  ln(Xi)  for  the  lognormal  distributions. 
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Table  1 .  Characteristics  of  Pores  on  I^fferent  Test  Planes  and  at 


jMfa|^Locj^onawhhmtibe225^T7451_PI|^Altej^ 


1  New  Pedigree  I 

Area% 

AR 

SI 

Q 

Pores/mnr^ 

No.  of  Fi^s 

N 

3.55 

0.114 

1.55 

0.109 

69.36 

234 

195 

3.67 

0.199 

1.49 

0.800 

0.478 

89.42 

215 

231 

1  ST.  center 

3.90 

0.136 

1.56 

0.805 

0.392 

69.36 

240 

200 

2.91 

0.071 

1.47 

0.831 

0.353 

64.34 

260 

201 

ILS.  surfiice 

2.73 

0.075 

1.52 

0.815 

0.650 

77.17 

220 

204 

3.07 

0.075 

1.53 

0.815 

0.546 

63.06 

260 

197 

3.91 

0.080 

1.52 

0.815 

0.163 

41.61 

200 

100 

3.50 

0.085 

1.52 

0.823 

0.465 

54.10 

200 

130 

4.03 

0.093 

1.46 

0.829 

0.270 

54.93 

200 

132 

1  Average 

3.47 

0.103 

1.51 

0.816 

0.381 

64.82 

225 

177 

1  OldPedij 

sree 

Bms] 

Ajrea% 

AR 

SI 

n 

% 

1 

No.  of  Fields 

N 

3.27 

0.125 

1.50 

0.811 

0.350 

95.71 

180 

207 

4.69 

0.368 

1.45 

0.834 

0.302 

158.68 

150 

286 

5.48 

0.286 

1.51 

0.820 

0.465 

92.94 

180 

201 

3.07 

0.071 

1.50 

0.813 

0.147 

65.33 

200 

157 

[EE5Z^3 

2.48 

0.058 

1.40 

0.845 

0.358 

76.04 

220 

201 

BSESilH 

2.73 

0.042 

1.41 

0.842 

0.238 

■Esm 

210 

128 

mmsm 

3.07 

0.142 

1.45 

0.834 

0.406 

130.60 

130 

204 

|LS,  quarter 

2.64 

0.042 

1.42 

0.839 

0.345 

50.77 

200 

122 

3.11 

0.080 

1.41 

0.843 

0.299 

68.66 

200 

165 

1  Average 

3.39 

0.135 

1.45 

0.831 

0.323 

%i.n 

186 

186 

LII-  Average  intercept  length. 
AR<  Aspect  ratio. 

SI  -  Sluq)e  ind«c. 

Q-  Orientation  factor. 

N  -  Numb«‘  of  detected  pores. 
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Table  2.  Parameters  of  Lognormal  Distnbution  Describing  Pore  Sizes  on 


^^^BdJPto>mji^^crtionswithin^tePtoeJ05^[74Slj\llo^ 


1  New  Pedigree  I 

Std.  Dev. 

R'^2 

N 

3.26 

1.55 

0.22 

0.9425 

195 

3.04 

1.74 

0.22 

0.9695 

231 

1  ST.  cento* 

3.65 

1.56 

0.22 

0.9450 

200 

3.04 

1.39 

0.19 

0.9374 

201 

LS.  sur&ce 

2.90 

1.38 

0.19 

0.9033 

204 

ST.  surfiwe 

2.23 

1.43 

0.20 

0.8601 

197 

3.62 

1.53 

0.30 

0.9293 

100 

3.53 

1.45 

0.25 

0.9492 

130 

SiFEi3 

3.18 

1.55 

0.26 

0.8980 

132 

1  Average 

3.16 

1.51 

0.23 

0.9260 

1  Old  Pedigree  I 

Std.  Dev. 

Avg.  Conf  95% 

R'^2 

N 

3.02 

1.45 

0.20 

0.9571 

207 

3.67 

1.74 

0.20 

0.9771 

286 

3.87 

1.63 

0.23 

0.8898 

201 

2.87 

1.37 

0.21 

0.8738 

157 

2.61 

1.36 

0.19 

0.9059 

201 

2.38 

1.42 

0.25 

0.8163 

128 

Eiill!S!!il 

2.69 

1.49 

0.20 

0.8895 

204 

2.54 

1.40 

0.25 

0.8369 

122 

1st.  quarto 

3.02 

1.43 

0.22 

0.8798 

165 

2.96 

1.48 

0.22 

0.8918 

N  -  Number  of  detected  pores. 
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Table  3.  Characteristics  of  the  True  3-D  Pore  and  C<»stituent  Particle  Sizes  fix  die  7050-T74S1  Plate  Alloys 


-  Largest  ellipsoid  radius. 

b.  a/c  -  Asjpect  ratios. 

e  •  Parameter  describing  goodness  of  fit. 

„  -  Average  size  on  indicated  planes  (fun). 


Table  4.  Characteristics  of  Consrituent  Particles  on  IMfferent  Test  Planes  and 


at  Diflferent  Locations  whbin  the  70S0~T74S1  Plate  Alloys 


1  New  Pedigree  I 

mmm 

Area% 

AR 

SI 

Cl 

Ptcls/mm^2 

No.  of  ndds 

N 

6.20 

0.542 

1.78 

0.748 

0.224 

94.32 

180 

204 

4.84 

0.480 

2.35 

0.663 

0.754 

69.29 

215 

179 

1  ST.  centor 

3.71 

0.845 

2.27 

0.663 

0.586 

297.77 

90 

322 

3.44 

1.125 

2.34 

0.674 

0.432 

509.76 

40 

245 

3.35 

0.745 

2.41 

0.652 

0.741 

371.19 

50 

223 

2.97 

0.458 

2.05 

0.712 

0.532 

295.45 

60 

213 

3.72 

0.350 

2.33 

0.671 

0.169 

140.79 

120 

203 

4.13 

0.638 

2.48 

0.659 

0.706 

216.39 

80 

208 

ST,  quarts 

3.77 

0.349 

^53 

0.646 

0.546 

135.72 

130 

212 

Avo^e 

4.01 

0.615  2.28 

0.676 

0.521 

107 

223 

i  Old  Pedi^ 

m _ 1 

Area% 

AR 

SI 

Cl 

Ptcls/mm^2 

No.  of  Fields 

N 

6.14 

0.629 

2.06 

0.704 

0.412 

94.32 

180 

204 

6.38 

0.524 

2.36 

0.655 

0.698 

66.91 

250 

201 

IP5ilt»!'J.T?4l 

4.57 

0.907 

1.62 

0.790 

0.560 

327.36 

60 

236 

4.04 

0.746 

1.71 

0.769 

0.318 

328.01 

51 

201 

3.64 

0.545 

1.90 

0.736 

0.632 

254.43 

70 

214 

3.12 

0.453 

1.76 

0.760 

0.348 

303.78 

60 

219 

EilBBES 

5.67 

0.839 

1.78 

0.761 

0.148 

186.80 

90 

202 

|LS,  quarter 

4.38 

0.708 

1.77 

0.761 

mmm 

249.68 

70 

210 

jKsm 

0.931 

1.87 

0.747 

252.06 

70 

212 

1  Average 

4.77 

0.698 

1.87 

0.743 

0.428 

229.26 

100 

211 

LII-  Average  intercept  length. 
AR-  Aspect  ratio. 

SI  -  Shape  index. 

Orientation  factor. 

N  -  Number  of  detected  particles. 
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Table  5.  Parameters  of  Lognormal  Distribution  Describing  Consdbient  Particle 
on  Indicated  Planes  and  Locations  within  the  Plate  70S0-T7451  Alloys _ 


1  New  Pedigree  I 

1  Sample 

Std.  Dev. 

Avg.Conf  95% 

R''2 

N 

5.30 

1.73 

0.24 

0.9885 

204 

ILS.  center 

5.10 

1.50 

0.22 

0.6633 

179 

3.45 

1.65 

0.18 

0.6627 

322 

3.76 

1.41 

0.18 

0.9151 

245 

3.78 

1.36 

0.18 

0.8964 

223 

3.41 

1.32 

0.18 

0.9076 

213 

4.13 

1.42 

0.20 

0.9215 

203 

4.29 

1.47 

0.20 

0.9408 

208 

4.07 

1.44 

0.19 

0.9408 

212 

4.14 

1.48 

0.20 

0.8707 

1  Old  Pedigree  I 

Std.  Dev. 

R'^2 

N 

5.27 

1.63 

0.22 

0.9611 

204 

5.80 

1.69 

0.23 

0.9926 

201 

3.81 

1.55 

0.20 

0.9588 

236 

4.07 

1.44 

0.20 

0.9556 

201 

4.01 

1.39 

0.19 

0.9641 

214 

3.34 

1.41 

0.19 

0.9312 

219 

4.74 

1.68 

0.23 

0.9428 

202 

il^EtlSSSI 

4.18 

1.52 

0.21 

0.9282 

210 

ST,  quarter 

4.62 

1.58 

0.9296 

212 

Average 

4.43 

1.54 

0.21 

0.9516 

N  -  Number  of  detected  particles. 
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Table  6.  Results  of  the  Tessellation  Analysis  for  Pores  and  Constituent 
Particles  in  the  Old  Pedigree  7050-T7451  Plate  Alloy 
OST  Plane,  Center  Section) 


] 

*ores 

BGS^i 

Minimum 

Maximum 

St  Dev. 

(jeom.Av. 

Median 

Nw-Ndglib. 

_ 

114 

2.83 

324 

65.3 

91.5 

116 

48.2 

2.83 

136 

33.2 

35 

42.5 

Local  V.F.  % 

.00302 

.000129 

0.136 

.0141 

.00168 

.00146 

Paiticie  Area  (lun^) 

30.8 

2.3 

320 

61.9 

14.1 

11.7 

Cdl  Area  (mm*) 

.0108 

.000439 

.033 

.00714 

.00842 

.00952 

CeU  Sides  No. 

5.74 

3 

15 

Particles 

Miniimim 

Maximiim 

St  Dev. 

GeontAv. 

Median 

Near-Ndgbb. 

_ 

38.0 

1.0 

194 

38.2 

44.0 

66.6 

23.2 

1.0 

119 

19.2 

17.9 

21.3 

jiES!!E!X3H 

.0413 

.000262 

4.38 

.292 

.007 

.00628 

Particle  Area  (um*) 

27.8 

3.00 

266 

44.5 

14.8 

12.2 

Cdl  Area  (mm*) 

.0033 

.000029 

.0034 

.0021 

.0022 

Oil  Sides  No. 

5.76 

3 

15 
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Table  7.  Values  of  Q  and  R  for  Pore  and  Particle  Distributions  in  7050-T7451  Plate  Alloy 
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N  -  Number  of  detected  features 
rdm  -  Random  distribution, 
new,  old  -  New  and  old  pedigree  alloys, 
obs  -  Experimental 


both  old  and  neTy^-TT^T^^f  different  plate  locations 
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Figure  3.  A  typical  pore  in  the  7050-T7451  plate  alloy. 


Figure  4.  A  collection  of  201  pores  found  on  the  ST-planes,  center  region,  old  7050- 
T745 1  plate  alloy.  Pores  were  collected  from  1 80  micrographs. 
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Figure  S.  The  average  pore  intercept  lengths  for  the  new  pedigree  70S0-T74S1  plate  alloy 
as  a  function  of  direction  at  the  surface,  quarts  depth  and  center  of  the  plate. 


Figure  6.  The  average  pore  intercept  lengths  for  the  old  pedigree  7050-T7451  plate  alloy 
as  a  function  of  direction  at  the  sur&ce,  quarter  depth  and  center  of  the  plate. 
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Figure  7.  Average  intercept  length,  Lg,  for  pores  and  particles  at  different  plate  locations 
fw  the  old  and  new  7050-T7451  plate  alloys. 
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Figure  8.  Cumulative  2-D  pore  size  distributions  for  the  70S0-T7451  plate  alloys. 
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Figure  9.  An  example  of  the  pore  size  distribution  on  LS  plane,  center  r^on,  old  7050- 
T7451  plate  alloy. 


Figure  10.  An  example  of  the  fit  of  the  lognormal  distribution  to  the  pore  size  data;  LS 
plane,  center  r^on,  old  70S0-T74S1  plate  alloy.  Solid  line  represents  lognormal 
distribution  and  open  squares  are  the  experimental  data. 
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N«w  pcdigna.  Quarttr 


Pora  Sin,  um 


New  ptdigrae,  turtle* 


Poi*  Sin,  lan 

Figure  11.  E}i|)eriinental,  2-D,  and  estimated  3-D  cunsilative  pore  size  distributions  in  the 
new  pedigree  70S0-T74S1  plate  alloy. 
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Figure  12.  ComfMuison  of  the  experimental  and  simulated  aumilative  2-D  pore  size 
distributions  in  the  new  pedigree  70S0-T74S1  plate  alloy. 


0  20  40  60  80  100  120  140  160  180 
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Figure  13.  The  angular  variation  of  the  normalized  average  intercqit  lengths  for  the  pores 
in  the  old  pedigree  7050-T7451  |^e  alloy,  on  ST  plane,  center  region.  See  text  for 
definitions  of  global  and  average  d's. 


Polar  Plot 
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Figure  14.  The  normalized  average  d  of  the  pores  plotted  in  polar  coordinates.  The  polar 
plot  gives  a  direct  visual  representation  of  the  pore  shapes. 
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Figure  IS.  Average  d»pe  indexes,  SI,  for  the  pores  and  particles  at  differoit  plate 
locations  for  both  old  and  new  pedigree  70S0-T74S1  plate  alloys. 
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Figure  16.  Average  orioitation  factors,  £2,  for  the  pores  and  particles  at  different  plate 
locations  for  both  old  and  new  pedigree  70S0>T74S1  plate  alloys. 


.  Figure  17.  A  collection  of 236  constituent  particles  found  on  ST  plane,  crater  region  of 
the  old  pettigree  70S0-T74S1  A1  plate  alloy.  Particles  were  collected  from  60  micrographs. 


Figure  18.  The  average  particle  intercept  lengths  for  the  new  pedigree  70S0-T74S 1  plate 
alloy  as  a  function  of  direction  at  the  sur&ce,  quarter  depth  a^  center  of  the  plate. 


Figure  19.  The  average  particle  intercept  lengths  for  the  old  pedigree  70S0-T74S1  plate 
alloy  as  a  function  of  direction  at  the  surfoce,  quarter  depth  and  coiter  of  the  plate. 
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Figure  20.  An  example  of  the  size  distribution  of  the  constituent  particles  on  ST  plane, 
colter  region  of  the  old  pedigree  70S0-T74S1  Al  plate  alloy. 


Figure  21.  An  ocample  of  the  fit  of  the  lognomal  distribution  to  the  particle  size  data;  ST 
plane,  coiter  region,  old  pedigree  7050-T74SI  Al  plate  alloy.  Solid  line  represents  the 
lognormal  distribution  and  open  squares  are  the  experimental  data. 
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Figure  22.  Cumulative  2-D  particle  size  distributions  for  the  7050-T7451  plate  alloy. 
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Figure  23.  The  average  variation  of  the  normalized  average  intercept  lengths  for  the 
constituent  particles  shown  in  Fig.  17.  See  text  for  defiiutions  of  global  and  average  d's. 


Polar  Plot 
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Figure  24.  The  normalized  average  d  for  the  constituent  particles  in  Fig.  17  plotted  in 
polar  coordinates.  The  polar  plot  gives  a  direct  visual  representation  of  the  particle  shapes. 
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Figure  25.  The  Dirichlet  tessellation  cells  constructed  on  the  pore  centers;  ST  plane, 
center  region,  old  pedigree  7050-T7451  plate  alloy.  Points  represent  locations  of  pore 
centers. 
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Figure  26.  The  normalized  average  intercept  lengths  of  the  tessellation  cells  constructed 
for  the  pores.  The  d-a  curves  for  a  rectanglar  and  random  point  patterns  are  shown  for 
comparison. 
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